Draft version March 22, 2013 

Preprint typeset using I^T^^X style emulatcapj v. 5/2/11 



o 



(N 



BIASES IN PHYSICAL PARAMETER ESTIMATES THROUGH DIFFERENTIAL LENSING MAGNIFICATION 
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ABSTRACT 

We study the lensing magnification effect on background galaxies. Differential magnification due 
to different magnifications of different source regions of a galaxy will change the lensed composite 
spectra. The derived properties of the background galaxies are therefore biased. For simplicity, we 
model galaxies as a superposition of an axis-symmetric bulge and a face-on disk in order to study the 
differential magnification effect on the composite spectra. We find that some properties derived from 
the spectra (e.g., velocity dispersion, star formation rate and metallicity) are modified. Depending 
on the relative positions of the source and the lens, the inferred results can be either over- or under- 
estimates of the true values. In general, for an extended source at strong lensing regions with high 
magnifications, the inferred physical parameters (e.g. metallicity) can be strongly biased. Therefore 
detailed lens modelling is necessary to obtain the true properties of the lensed galaxies. 
Subject headings: Gravitational lensing; high redshift galaxy 



o 

u 

Oh! 

o 



cn 
(N 

rn 
O 
m 



1. INTRODUCTION 

The properties of high redshift objects, mostly galax- 
ies, provide important information on the early evolution 
of galaxies and the role they played in the cosmic reion- 
ization (e.g. Loeb 2010; Fan 2012). Observations of the 
high redshift galaxies are challenging due to large dis- 
tances and faint magnitudes. With the growing power of 
telescopes, the number of high redshift galaxies found is 
increasing rapidly, either by deep imaging or the color se- 
lection (the y-band, r-band "dropout" technique) (e.g. 
Stanway et al. 2004; Hildebrandt et al. 2009). More- 
over, the F— band and J— band dropouts can potentially 
detect objects at redshift 2 ~ 8 (e.g. Yan et al. 2012). 
Galaxy clusters, as nature telescopes can enhance the 
ability to detect high redshift galaxies (Soucail 1990). 
The method has been implemented to study galaxies over 
a wide range of redshifts (e.g. Ellis et al. 2001; Bradac 
et al. 2009; Hall et al. 2012) as well as Lya spheres (Li 
et al. 2007). It has been shown that the search for the 
high-z galaxies (z > 7) in a lensing field is more efficient 
than that in blank fields, and the maximum efficiency is 
reached for lens clusters at 2; ~ 0.1 — 0.3 (Maizy et al. 
2010). 

Gravitational lensing provides a direct way for study- 
ing the mass distribution of the large scale structures 
in the universe as well as galactic- and cluster-sized ha- 
los (e.g. Bartelmann & Schneider 2001; Treu 2010). On 
the other hand, the light from distant galaxies can be 
magnified by several orders of magnitude by the gravi- 
tational potential well of massive galaxy clusters. The 
effective solid angle of the survey volume decreases in 
the same way. However, since the luminosity function is 
exponential at the bright end, the magnification signifi- 
cantly increases the number counts of luminous galaxies 
(e.g. Lima et al. 2010; Er et al. 2013). In addition, the 
lensing magnification improves the spatial resolution of 
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distant galaxies (e.g. Stark et al. 2008; McKean et al. 
2011; Fan et al. 2012; Jones et al. 2013). For muffiply- 
imagcd systems, the positions and magnitudes of the im- 
ages depend on several factors, e.g., lens and source red- 
shifts, lens mass profiles (including substructures, e.g. 
Mao & Schneider 1998) and intrinsic properties of the 
background sources. Thus gravitational lensing provides 
us with more information for studying both lens galaxies 
and background sources. 

Lensing magnification itself is achromatic. However, 
magnification varies with the position of the source. The 
spatial profiles at different wavelengths may not be the 
same. Thus a source at different wavelengths may be 
magnified differently (Blandford & Narayan 1992). It 
has been noticed that differential magnification can bias 
the derived results using the spectral energy distribution 
(SED) of lensed sources (e.g. Blain 1999; Pontoppidan & 
Wiklind 2001), or using the ratio of spectral lines (e.g. 
Downes et al. 1995). Hezaveh et al. (2012) modelled their 
sources with two components (a compact and a diffuse 
one) and showed that the compact one will usually be 
magnified by a larger factor than the diffuse one in strong 
lensing. 

In this paper, we will study the bias in the derived 
spectral physical parameters due to differential magnifi- 
cation. We model the background galaxy as a sum of 
a bulge and a disk with different sizes and spectra. Us- 
ing ray-tracing simulations, the lensing magnifications of 
different regions are calculated and applied to obtain the 
composite observed spectrum. We start with a discussion 
of the basic formalism in section 2, present our model and 
magnification effects in section 3, and then discuss our 
results in section 4. The cosmology adopted here is a 
ACDM model with parameters based on Kilbinger et al. 
(2012): f^A = 0.718, f7„i = 0.286, erg = 0.804, a Hubble 
constant = lOO/i kms^^ Mpc^^ with h = 0.693. 

2. LENSING MAGNIFICATION 

The fundamentals of gravitational lensing can be found 
in Bartelmann & Schneider (2001). The thin- lens ap- 
proximation is adopted in this paper, implying that the 
lens mass distribution can be projected onto the lens 
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Fig. 1. — The simulated image of a background source (left) and 
the lensed (right) images. Notice that the scales and units of the 
two panels are different. In the right panel, the lens is at the origin 
of the map with 6e = 7.8 arcsec (~ 15kpc). The solid line is the 
critical curve of the lens. The lens and source are at redshift 0.1 
and 2.0 respectively. 
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Fig. 2. — The left panel shows the surface intensity profile of 
the source galaxy: an axis-symmetric bulge, a face-on disk compo- 
nent and the sum. The right panel shows the velocity dispersion 
variation as a function of radius. 

plane perpendicular to the line of sight. We denote an- 
gular coordinates on the lens plane as 0, and those on 
the source plane as (3. The lens equation can be written 
as 

I3 = e-^a{e) = 0-cy{e), (1) 

where Z?ds and Ds arc the angular diameter distances 
from the lens to the source and from the observer to 
the source respectively. The deflection angle a can be 
calculated from the lens model. 

We denote the brightness distribution of the source by 
/"(/S). Since lensing conserves the surface brightness, 
the brightness distribution of the image is thus I{d) = 
P{(3{6)). The total flux of the source and image are 
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d2/3r(/3), 



S 



(2) 



respectively. The magnification is defined as fi = S/Sq. 

3. LENSING EFFECTS ON THE SPECTROSCOPIC 
PROPERTIES OF BACKGROUND SOURCES 

Lensing magnification varies with the relative positions 
of the lens and source, especially when the source is 
close to the caustics. For extended galaxies, different 
components differ in intrinsic sizes, thus their magnifi- 
cations will be different. Therefore the observed spec- 
trum of the background galaxy is changed by lensing. 
The estimation of galaxy properties differs from its in- 
trinsic ones, like velocity dispersion (cr*), star formation 
rate (SFR) and metallicity of galaxies (Z). In Hezaveh 
et al. (2012), the galaxy is modelled as the sum of an 



extended component and an inner core. The total flux 
is the sum of the two components and the total mag- 
nification is their weighted average. A similar approach 
will be adopted here. In reality, galaxies are structurally 
complicated, e.g. they host some clumpy star formation 
regions. Moreover, their morphology may not be regu- 
lar, especially during the early stage of galaxy formation. 
The simple model here is for illustration purposes. 

3.1. Light profile of the background galaxy 

The luminosity profiles of the bulge and the disk are 
described as follows: 



/(i?)=/eexp{-7.67 {R/ReY^^-l } 
I{R)^IdCxp{-R/Rd), 



(3) 
(4) 



where /g is the surface brightness at the bulge effective 
radius Re, Id is the central surface brightness and Rd is 
the scale length of the exponential disk. The bulge to 
total ratio is (Binney & Tremaine 1987) 
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We use Id ~ 1.0 and Ig ~ 0.44 (the unit is arbitrary 
here, since we only need the luminosity ratio between 
the bulge and the disk). With a small bulge assump- 
tion in this study (B/T = 0.2) (Lackner & Gunn 2012), 
we take Re = 0.7kpc and Rd = 1.7kpc, corresponding 
to 0.08 and 0.2 arcsec at redshift Zg ~ 2 respectively. 
For simplicity, we assume an axis-symmetric bulge and 
a face-on disk (see the left panel of Fig. 1). 

For convenience, we will use the angular separation 
in the rest of this paper (1 arcsec equals to 8.7 kpc in 
the source plane Zg = 2, and 1.9 kpc in the lens plane 
Zd = 0.1). The velocity dispersion is estimated from the 
surface brightness using the isothermal sheet. Further 
approximation, assuming a uniform disc scale height, the 
velocity dispersion of the source as a function of radii 
is given by ct»(6') oc \/I{6) (Binney & Tremaine 1987, 
Chapter 4.4), and we set the central value as 218kms~^ 
(see the right panel of Fig. 2). 

The stellar velocity dispersion and the metallicity of 
source galaxy relies on stellar continuum, while the star 
formation rate (SFR) can be estimated from the [O 
ii] emission lines. Hence we simulate the spectra of bulge 
and disk by combining the emission-line spectra and the 
stellar continuum. For stellar continuum, four artifi- 
cial stellar populations are used to generate the spectra 
(Bruzual & Chariot 2003). The weight of each stellar 
population is set to be equal (1/4). Due to the metallic- 
ity gradient of galaxies (e.g. RoUeston et al. 2000), the 
metallicity of the bulge is set to be larger than that of 
the disk region (see Table 1 for more details on the stellar 
populations of the bulge and disk) . 

The emission-line spectra of the bulge and disk are 
also simulated separately. For the bulge there are two 
components: a narrow-line region (NLR) of the central 
active galactic nucleus (AGN) and a starburst in the 
bulge. For the disk there is only a starburst compo- 
nent. In order to study the stellar continuum, we can 
only take into account the emission of a type II AGN 
by assuming that the emissions from the accretion disk 
and the broad line region are obscured by the torus. 
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TABLE 1 

Stellar populations used for the bulge and disk 



Name 


Age 


Mctallicity ZIZq 




Bulge 




agcll6_m42 


1.0151!E+08 


0.2 


age070_m62 


l.OOOOE+07 


1.0 1 


agel35_m62 


9.0479E+08 


1.0 " 


agel50.m72 


2.5000E+09 


2.25 


Disk 


agc055_m42 


5.0100E+06 


0.2 


age070.m42 


l.OOOOE+07 


0.2 


agel39.m42 


1.4340E+09 


0.2 


agell6_m62 


1.0152E+08 


1.0 



The type II AGN emission-line spectrum is taken from 
Reyes et al. (2008). The starburst emission- Une spectra 
are taken from SDSS DR7 with the selection criteria of 
[OIII]/i/^ - 0.35 and [NII]/ffa 1.0, which is appro- 
priate for the star-forming gas rich galaxies at redshift 
2. The SFRs of the bulge and disk, which are estimated 
by using [O ii] luminosities, are set to be proportional to 
the surface brightness. 

With the chosen spectra and luminosity profiles of the 
bulge and the disk, the lensed spectrum of the back- 
ground galaxy can be obtained after knowing the mean 
magnification of each part of the source: 



F(v, e) = tibFb{v, 6) + iidFd{y, 9). 




(6) 



where Hb and Hd are the magnifications of the bulge re- 
gion (< 0.1 arcsec) and the disk region, which can be 
calculated from the lensing model. 

3.2. Lens modelling 

In order to study the differential magnification effect 
on the spectrum of background sources, we perform ray- 
tracing simulations. For simplicity, we model our lens 
as a singular isothermal ellipsoid with a velocity disper- 
sion of 770 km/s (Af200 ~ 1 x IO^'^/t.^^Mq) and cUipticity 
e = 0.2. The lens is placed at redshift Zd = 0.1, which 
gives an Einstein radius 6'e — 7.8 arcsec. We map a grid 
of pixels from the lens plane to the source plane using 
the lens equation to obtain the surface brightness distri- 
bution of the image. From this distribution, the total 
flux of the image is obtained using Eq. 2. We will cal- 
culate the magnifications separately for the inner region, 
(< 0.1 arcsec), and for the outer region (/id)- 

The magnification ratio between the bulge and the disk 
strongly depends on the relative position of the source 
and lens, and the intrinsic size of source galaxies. In the 
left panel of Fig. 3, one can see that when the bulge of 
the source galaxy is close to the caustics, the magnifica- 
tion differences in the bulge and the disk are significantly 
larger. An interesting point is when the disk is closer to 
the caustics than the bulge, the ratio becomes smaller 
than 1 (/id > /if,). When the bulge crosses the caus- 
tics, the reverse occurs {^b > l^d)- As a consequence, a 
slight shift of the source position (or a change of the lens 
model) may cause a significantly different magnification 
ratio. One can also see that the region where the differ- 
ential magnification is strong, the absolute magnification 
is also high (e.g. /x > 5, right panel of Fig. 3). 

We first place the source at (3 = (0, 3) arcsec on the 
source plane (and study another position (0.3, 2.5) in the 
end) . Two images are generated by lensing (see the right 
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Fig. 3. — Left : The ratio of fih/Hd on the source plane, given 
on linear grey scale and contours. The blue dashed lines enclose 
the region where fit/fJ-d < 0.85. The green lines enclose the dark 
region where the ratio is greater than 1. The red solid line is the 
caustics. The cross and plus represent the positions of our mock 
source galaxies. Right : The mean magnification of the whole 
galaxy on the source plane. The green dashed line encloses the 
region where |^| > 5. 
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Fig. 4. — The luminosity- weighted mean velocity dispersion as a 
function of the fiber size. 

panel of Fig. 1). We only study the primary image in this 
paper, since it is more luminous and easier to detect. The 
mean magnification of the bulge (disk) region is /if, = 15 
(/id = 10). 

We simulate the velocity dispersion cr* detected by dif- 
ferent sizes of fibers. The centre of the fiber is aligned 
with the centre of the lensed source. The mean veloc- 
ity dispersion is weighted by luminosity. In Fig. 4, the 
solid and dashed lines represent the a^,{9) after and be- 
fore lensing respectively. When the fiber size becomes 
very small, the measured velocity dispersion approaches 
the intrinsic ones; in the other cases, the measurements 
of the lensed galaxy are larger than the initial ones for 
/3 = (0,3). 

3.3. Results from spectral analysis 

The (T, , SFR and metallicity (Z) variation of the source 
and the lensed galaxy can be obtained through the anal- 
ysis of the composite spectra. According to the source 
model and the parameters listed in Table 2, we can fit 
the galactic continuum using STARLIGHT (Cid Fernan- 
des et al. 2005). This code employs 45 stellar population 
templates which are composed of 3 metallicities and 15 
ages (Bruzual & Chariot 2003). The differences between 
the spectra of the lensed image and the initial galaxy are 
listed in Table 2. The velocity dispersion of the lensed 
image is higher in the case of jib/Hd = 15/10, since the 
velocity dispersion of the bulge is larger than the disk. In 
order to get a reliable result, we fit the spectra by using 
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the direct pixel fitting method (Greene & Ho 2006; Ho 
et al. 2009; Ge et al. 2012), and obtain the same results. 

With the continuum subtracted, we can obtain the 
emission-line only spectra. The SFR is estimated by 
fitting [OHJA3727. Here the SFR is normahzcd to 
lOMoyr-i (e.g. Santini et al. 2009). The SFR derived 
from the lensed image is enhanced by about 11 times in 
the case of fib/fJ'd = 15/10. After correcting by the mean 
magnification of the whole galaxy (/i = 12.7), the SFR 
result is smaller than the true value by ^ 10%. 

The metallicity of the whole galaxy is again obtained 
from the STARLIGHT, Z = Y,Xj x Zj, where x-j is the 
weight of the j-th stellar population, and Zj is its metal- 
licity. The result from the lensed image is significantly 
larger than the initial one. The reason is that both the 
metallicity and the magnification of the bulge arc higher 
than those of the disk, and the weight of the bulge is 
enlarged by lensing. 

To see the sensitivity of the results to the source posi- 
tion, we place the source at a different location (0.3, 2.5). 
The results are given in the column 'lensed 2' in Table 2. 
The magnification ratio is different from the first case: 
fJ-b/ fJ-d = 7/10. In this case, the contribution of the bulge 
to the total lensed spectrum is suppressed, which leads 
to a smaller cth. and metallicity while the SFR is overes- 
timated. 

When we observe the lensed galaxy with fibers, the 
light may actually come from different parts of the bulge 
and the disk, which makes the measurement of velocity 
dispersion, SFR and metallicity of the source more com- 
plicated. We perform simulations with a fiber size of 3 
arcscc diameter. The centre of the lensed image is iden- 
tified as the brightest position of the image. The results 
are given in the bottom part of Table 2. For both source 
positions, similar biases are obtained. 

4. SUMMARY AND DISCUSSION 

In this paper, we have studied the differential magni- 
fication effect on the spectrum of a lensed background 
galaxy. Two components (a bulge and a disk) with dif- 
ferent luminosity profiles and spectra are used to model 
the background galaxy. Ray-tracing simulations are em- 
ployed to calculate the lensing magnification. We find 
that the derived properties of the lensed galaxy are 
changed, e.g. SFR, metallicity and velocity dispersion. 
Velocity dispersion and metallicity can either increase or 
decrease, depending on the relative position of the source 
and the lens. The SFR is strongly enhanced but can be 
corrected by dividing the mean magnification. However, 
a slight bias of SFR still remains after correction. Not 



surprisingly, the differential magnification will be strong 
when the source is close to the caustics, and the absolute 
value of the magnification is also high (e.g. \fj,\ > 5). In 
this case, a slight shift of the source position will cause 
a dramatic change in the derived results. Therefore, one 
needs to take into account the differential magnification 
effect for highly magnified extended sources. 

Because of observational limitations, the black hole vs. 
velocity dispersion (A/bh — cr*) relation at high redshift 
(e.g. Tremaine et al. 2002) is uncertain. Gravitational 
lensing boosts our chances of studying black holes at high 
redshift. However, we found that the measured ct, may 
be an unreliable indicator of the black hole mass. In 
practice, for very high-redshift galaxies, it is easier to 
observe with slits than fibres since the background sub- 
traction is more reliable. To avoid extra bias entering 
the results, one needs to consider carefully the region of 
the background galaxy covered when taking a spectrum. 

Our study is simplistic in several ways. For example, 
source galaxies may be irregular, lacking regular bulge 
and disk components. Additionally, high redshift galax- 
ies can be complicated, e.g. they may have clumpy star 
forming regions with different spectra (Genzel et al. 2008; 
Jones et al. 2013). This will mislead the calculation of de- 
rived physical properties. Moreover, lens clusters arc also 
complicated. The magnification map is strongly affected 
by the shape and mass profile (including substructures) 
of the lens. Therefore accurate lens modelling is neces- 
sary when deriving the properties of lensed galaxies. 

On the other hand, multiple images can provide more 
constraints on both the lens and source. With Integral 
Field Unit (IFU) spectra for different components of the 
multiple images at high redshift, we can study the lens 
mass distribution and the intrinsic properties of back- 
ground galaxies in greater detail. With iterative mod- 
elling of the lens and source, a relatively precise source 
image can be reconstructed. The bias due to differential 
magnification can be strongly suppressed. 

Current telescopes, such as the Hubble Space Telescope 
(HST) or the Keck telescope can already study differen- 
tial magnifications. The next-generation telescopes with 
adaptive optics, such as the Thirty Meter Telescope and 
the E-ELT, will allow even more accurate studies of dy- 
namics in high redshift galaxies. 
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TABLE 2 

Results of differential magnification 



Parameters 


Bulge 


Disk 


lenscd 1 


lenscd 2 


Source 


total flux 




15(7) 


10(10) 


15+10 


7+10 


1+1 


luminosity 


1 


4 


58 


50 


5 


cr« (km/s) 


131 


48 


67 


58 


59 


SFR (10/i-lMo/yr) 


2 


8 


113(8.9) 


93(11) 


10 


metallicity Z/Zq 


0.95 


0.4 


0.59 


0.47 


0.51 


flux within 3 arcscc diameter flbcr 




15(6.7) 


13.5(6.6) 


15+13.5 


6.7+6.6 


1+1 


luminosity 


1 


0.67 


19 


11 


1.67 


(T» (km/s) 


164(134) 


53(52) 


118 


76 


110 (74) 


SFR (lOh^iMo/yr) 


6 


4 


144(9.6) 


67(10) 


10 


metallicity Z/Zq 


0.95 


0.4 


0.8 


0.58 


0.61 



Note. — The column 'lensed 1' represents the lensed results for the source position (0,3) (black cross in Fig. 3), while the column 
'lensed 2' represents the results for the source position (0.3, 2.5) (purple plus in Fig. 3). cr* is the velocity dispersion estimated using the 
method in Gc et al. (2012). The SFR in brackets is corrected by the mean magniflcation of the whole lensed galaxy. 
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